Soliton-sound interactions in quasi-one-dimensional Bose-Einstein condensates 
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Longitudinal confinement of dark solitons in quasi-one-dimensional Bose-Einstein condensates 
leads to sound emission and reabsorption. We perform quantitative studies of the dynamics of a 
soliton oscillating in a tight dimple trap, embedded in a weaker harmonic trap. The dimple depth 
provides a sensitive handle to control the soliton-sound interaction. In the limit of no reabsorption, 
the power radiated is found to be proportional to the soliton acceleration squared. An experiment 
is proposed to detect sound emission as a change in amplitude and frequency of soliton oscillations. 
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PACS numbers: 03.75.Lm, 42.65.Tg 

The experimental realisation of dilute atomic Bose- 
Einstein condensates (BEC) pj has introduced an un- 
paralleled platform from which to study the dynamical 
behaviour of nonlinear systems. Of great interest is the 
stability of topological structures, such as vortices and 
dark solitary waves. Recent experiments at low tem- 
peratures have highlighted the temperature-independent 
crystallization and decay Q of vortex structures. Re- 
lated theoretical work has indicated the importance of 
vortex-sound interactions Q. Following the recent ex- 
perimental observation of dark solitary waves (henceforth 
referred to as dark solitons) in BEC's H, one can now 
study the coupling of solitons to the background sound 
field. 

Dark solitons are robust localised one-dimensional de- 
fects characterised by a notch in the ambient condensate 
density and a phase slip across the centre. They are sup- 
ported in repulsive (defocussing) nonlinear media where 
the kinetic dispersion of the wave is balanced by the non- 
linear interaction. In three-dimensional geometries they 
are prone to decay into lower energy topological struc- 
tures such as vortex rings [(j. However, in quasi-one- 
dimensional (quasi-lD) geometries, where the transverse 
condensate size is of the order of the healing length, soli- 
tons are expected to exhibit longer lifetimes, due to the 
suppression of transverse excitations . 

In the context of nonlinear optics @, dark solitons 
on a homogeneous background can become intrinsically 
unstable to modifications of the nonlinearity within the 
optical medium (e.g. due to saturation effects) [j|. In 
quasi-lD BEC's, modified nonlinearities can arise from 
the transverse dimensions that have been integrated out, 
as considered for longitudinally homogeneous conden- 
sates in [To) . The longitudinal confinement of current 
atomic BEC experiments leads to an additional mech- 
anism by which the integrability of the system can be 
broken [n], 0, ^| . In this paper we perform detailed 
quantitative analysis of the stability of dark solitons in 
quasi- ID geometries under longitudinal confinement, and 
find that dissipation due to the longitudinal inhomogene- 
ity dominates over decay arising from coupling to trans- 
verse modes [T^. Our axial configuration consists of a 
tight dimple trap (where the soliton resides) embedded 
in a weaker harmonic potential. Such a geometry en- 



ables us to control and examine the interplay between 
sound emission and reabsorption by the oscillating soli- 
ton. In the limit of no reabsorption, the rate of emission 
is found to be proportional to the local soliton acceler- 
ation squared, as predicted for optical solitons in homo- 
geneous nonlinear waveguides @. Finally we propose a 
controllable experiment where the soliton-sound interac- 
tions can be quantified by measuring the amplitude and 
frequency of the soliton oscillations. 

Our analysis is based on numerical simulations of the 
Gross-Pitaevskii Equation (CPE), describing the dynam- 
ics of weakly-interacting BEC's in the limit of low tem- 
perature, 
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Here ip is the order parameter of the system, V ex t the con- 
fining potential, m the atomic mass and g = AmY?ajm 
the scattering amplitude, where a is the s-wave scat- 
tering length. In the limit of tight transverse confine- 
ment, the condensate dynamics are effectively longitudi- 
nal [Tj|. By integrating out the transverse degrees of 
freedom, the GPE reduces to a one-dimensional (ID) 
form with a modified coupling coefficient g 1D — g/(2irl 2 ± ), 

where l ± = \J h/mui ± is the transverse harmonic oscilla- 
tor length and lo ± the transverse trapping frequency. On 
a homogeneous background density n 1D , a dark soliton 
of speed v and position, z' = (z — vt), is defined as, 

tp(z,t) = v^[[3t&nh(f3z'/0 + i(v/c)}e-^» t/h (2) 

where /i 1D = g 1D n 1D is the one-dimensional chemical po- 
tential, c = yj /! 1d /to the Bogoliubov speed of sound, 
[3 = y/l — (v/c) 2 , and £ = ft/^/m/x 1D is the healing 
length characterising the size of the soliton. 
We consider a trap of the form 
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where we have defined the 'cut-off' depth of the inner 
'dimple' trap as Vo — \mJ\z 1 (see Fig. 1). Sound waves 
generated by the oscillating soliton can either be trapped 
within the dimple or escape to the outer trap, depend- 
ing on the value of the cut-off Vq relative to the chemical 
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FIG. 1: (a) Initial longitudinal condensate density (solid line) 
for a double trap geometry (dashed line) with a soliton centred 
at the origin, (b) Enlarged image of central region. 
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potential /i 1D . For a sufficiently weak outer trap, the 
density in the outer region and close to the dimple be- 
comes approximately homogeneous (Fig. 1(b)), and as a 
first approximation we initially solve the ID GPE in the 
limit uj^ = 0. We shall return to the experimentally rele- 
vant case of a>£ 7^ and 3D simulations later. The effect 
of the emitted sound on soliton dynamics can be deter- 
mined by monitoring the 'soliton energy' a quantity 
determined by integrating the functional, 
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across the soliton region (z s ± 5£) and subtracting the 
corresponding contribution of the background fluid. 

Fig. 2 shows the temporal evolution of the soliton 
energy for different cut-offs (corresponding to different 
amounts of sound reabsorption), all of which display a 
periodicity of approximately u z l due to the motion 
of the soliton in the trap [ll|, [12, uM ■ One can discrimi- 
nate two limiting cases: (i) For sufficiently high cut-offs 
(here Vq > 1.2/i 1D ) all emitted sound becomes trapped 
in the dimple, leading principally to the excitation of the 
dipole mode [ll| • The combined effect of the dipole mode 
and soliton, oscillating at different frequencies, induces 
beating of characteristic period r w 22 (see Fig. 2). (ii) 
For low cut-offs (Vq < 0.4/i 1D ) there is essentially pure 
emission (no reabsorption) of sound by the soliton. This 
behaviour is consistent with the change in energy experi- 
enced by a classical particle of negative mass oscillating 
in a harmonic trap, subject to a dissipative force which 
induces an energy loss proportional to the acceleration 
squared 0|. The observation that the soliton dynam- 
ics changes significantly around Vq ~ /x 1D (since sound 
excitations have energy of the order of the chemical po- 
tential) provides a sensitive handle on the soliton-sound 
interaction, which can be controlled experimentally. 

The dissipative dynamics of a dark soliton in a shal- 
low dimple trap are investigated in Fig. 3. As the soli- 
ton oscillates, it becomes asymmetrically deformed [Tfil ] 
and tries to adjust to the inhomogeneous background by 
radiating counter-propagating sound pulses (Fig. 3(a)). 
Fig. 3(b) shows the power emitted by the soliton as a 
function of time. At the bottom of the trap there is 
no sound emission since Eq. (2) becomes the local so- 



FIG. 2: Soliton energy versus time for a soliton of initial speed 
0.5c and position z = in a relatively weak inner harmonic 
trap (/i 1D = 70hui z ) for various potential cutoffs Vb (from top 
to bottom): l.l/i 1D , 1.02/i 1D , l/i 1D , 0.8/i 1D and 0.4^i 1D . The 
lowermost dotted line corresponds to the theoretical predic- 
tion of Eq. (5), valid in the limit of no reabsorption. 



lution to the GPE, while at the zenith of the first few 
oscillations, the soliton radiates maximum power. This 
sound emission (which modifies the soliton depth) leads 
to an increase in the soliton speed and amplitude of os- 
cillations, which in turn causes a gradual increase in the 
peak power radiated. As the soliton approaches the edge 
of the dimple, it experiences a smoothing of the effective 
potential due to fluid healing. This reduces the emission 
at the trap edge, as manifested in the 'intermediate' dips 
at u z t ~ 25 and 29.5. Eventually at uj z t ss 35 the soliton 
escapes the trap [Tl| . 

In the context of nonlinear optics, multiscale asymp- 
totic techniques predict that an unstable optical soli- 
ton, propagating in a homogeneous nonlinear waveguide, 
emits sound at a rate proportional to the soliton acceler- 
ation squared 0], such that, 
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Here n is the local background density, S s the total phase 
slip across the moving soliton, and iV s = J (n — IV* | 2 ) dz 
the number of particles displaced by the soliton. In the 
optical case, the breakdown of integrability leading to 
sound emission arises from modifications in the nonlin- 
earity within the medium, as opposed to the longitudinal 
confinement discussed above. It may therefore come as a 
surprise that this expression could also be applicable to 
the inhomogeneous case. In fact, this can only be valid 
in double trap structures of Eq. (3) if the dimple trap 
is sufficiently shallow that it essentially allows complete 
escape of sound to the outer trap; even then, this re- 
sult will only hold until the emitted sound returns to the 
dimple region after reflection from the outer trap. In this 
limit, the asymptotic conditions become practically indis- 
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FIG. 3: (a) Space-time plot of a soliton oscillating in a dim- 
ple trap with Vo = 0.4/^ 1D where fi 1D = 35hu z . Deviations 
from background density correspond to sound (light/dark re- 
gions indicate high/low densities), with the sound amplitude 
roughly 2% of the background, (b) Power emitted by the 
soliton as calculated from the energy functional (solid line), 
compared with the prediction of Eq. (5) (circles) . 




FIG. 4: (a) Soliton energies obtained from the 3D GPE (thin 
lines) vs. corresponding ID GPE results (thick lines) for cut- 
offs of Vo = 1.1/j, 1d (top) and Vo = 0Afi 1D (bottom). The 
longitudinal parameters are as in Fig. 2, and trap aspect ratio 
is set to lj±/uJz w 250, such that the initial density profile 
in the ^-direction is identical to the one-dimensional density, 
(b) The oscillation in the total longitudinal energy due to the 
coupling to transverse modes, as measured by integrating the 
energy functional longitudinally. 



tinguishable from the homogeneous ones, thus justifying 
the above approach. The validity of Eq. (5) is not depen- 
dent on the origin of the instability, but rather on the fact 
that the emitted sound escapes from the soliton region. 
Note that the effect of the dimple confinement is implicit 
in the variations of the parameters S s , c and n, and, con- 
trary to the homogeneous case 9] , this result will hold for 
all soliton speeds, due to the continuously-induced 'soli- 
ton instability' caused by the longitudinal confinement. 
Fig. 3(b) confirms the applicability of Eqs. (5) (circles) in 
the limit Vq <C /x 1D by direct comparison to the ID GPE 
simulations (solid line). Further evidence is provided in 
Fig. 2 which shows clearly that as the cut-off is decreased, 
the variation of energy with time (solid lines) approaches 
the above analytical result in the limit of no reabsorption 
(dotted line). We stress that Eq. (5) will not hold in the 
usual infinite harmonic traps, where multiscale asymp- 
totic technique must be combined with boundary layer 
theory 

Simulations of the 3D GPE with cylindrical symme- 
try (still under the assumption of a homogeneous outer 
trap) confirm the qualitative findings of Fig. 2. In partic- 
ular, in the limit /j, 3D /huj ± sa 1.2 (where the soliton was 
recently predicted to be stable against transverse decay 
|lf)j | ) we find the dominant decay mechanism to be due to 
axial inhomogeneity. The coupling to transverse modes 
|13| merely leads to an additional oscillation of the soli- 
ton energy in time (Fig. 4), whose amplitude increases 
as the transverse confinement is relaxed 0|- Subtract- 
ing the average oscillation in the 3D longitudinal energy 
(Fig. 4(b)) from the energy of the soliton region gives es- 



sentially perfect agreement with the ID results (Fig. 2). 

Finally we extend our earlier treatment to the experi- 
mentally relevant case of uj^ ^ 0. We consider a BEC in a 
quasi-lD geometry |l8j in a cylindrically symmetric trap 
with frequencies uj^ (longitudinal) and lo±_ (transverse), 
and then add an inner dimple potential (uj z ), such that 
the longitudinal confinement is as shown in Fig. 1. This 
could be realized, for example, by a dimple microtrap 
created optically within a larger magnetic trap. The key 
points for such a realization are: (i) a very weak longi- 
tudinal trap frequency u)^, such that the sound escaping 
from the inner trap is not reflected back from the outer 
trap within the timescale of interest, and (ii) a sufficiently 
tight inner trap lu z such that the oscillating soliton emits 
sound at a rate fast enough that its amplitude increases 
noticeably before other decay mechanisms become im- 
portant. 

Experimentally, the effect of soliton-sound interactions 
can be probed by determining the soliton path via re- 
peated time of flight measurements, and comparing paths 
for different depths of the dimple trap. The best signa- 
ture for this effect arises when comparing the case of an 
effectively infinite cut-off to that of a lower cut-off, which 
nonetheless permits a sufficient number of soliton oscil- 
lations. Since the soliton rapidly escapes the trap for 
Vo < /U 1D , we thus propose the use of Vq = fi 1D as a low 
cut-off, for which the soliton remains confined and ulti- 
mately decays to a sound pulse. Such a comparison has 
been performed for the double-trap system in Fig. 5, for 
both harmonic and gaussian dimple traps. One clearly 
observes a change in both amplitude and frequency of the 
soliton oscillations as a result of sound emission. These 
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FIG. 5: Soliton path in a longitudinal double-trap configura- 
tion for (a) tight harmonic inner trap and (b) gaussian dimple 
Vo[l — exp(— az 2 )]. Each graph shows the cases of a cut-off 
= Mm (dotted line) compared to an effectively infinite 
cut-off Vo = 5^i 1D (solid line), for a soliton with initial speed 
0.2c. The gaussian dimple in (b) is chosen such that it yields 
same initial oscillation frequencies as the harmonic dimple 
(achieved for a — u 2 /2Vo). The outer trap is harmonic in 
both cases, with relatively weak confinement ljq — u z /10. 



effects become more noticeable in the gaussian trap, for 



which the required cut-off for total sound entrapment is 
larger than the harmonic one. Choosing 'typical' param- 
eters ujq — 2tt x 5 Hz, lu z = lOw^ and uj ± — 250u>q and 
M.-SD = 8^57 where Zo = (ujqui 2 ^) 1 ' 3 , we find respectively for 
a 23 Na ( 87 Rb) condensate, a one-dimensional peak dim- 
ple density of n 1D = 5 x 10 7 (1.5 x 10 7 ) m _1 , condensate 
atom number TV = 18, 000 (3,500) and a total observa- 
tion timescale for Fig. 5 of r exp = 220 ms, consistent with 
the expected soliton lifetime in such systems 0, [n|, 112 ■ 

In summary, we have performed quantitative studies 
of the dynamics of dark solitons oscillating in trapped 
quasi-one-dimensional Bose condensates under longitudi- 
nal double-trap confinement. This geometry is ideal for 
controlling the interplay between sound emission and re- 
absorption, and we have identified two limiting cases: (i) 
In the limit of no reabsorption, the power emitted by the 
soliton is proportional to the square of its acceleration, 
(ii) The opposite limit of effectively infinite traps (such as 
harmonic traps in which soliton experiments have been 
performed) leads to the anticipated stabilization against 
soliton decay. We suggest that the sound radiation can 
be quantified by tracking the position of the dark soli- 
ton in suitably engineered traps, which allow for precise 
control of the soliton-sound interaction. 
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